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Abstract 

We present a time-dependent analysis of CP violation in p'^vr^ decays based on a 140 fb~"^ 

data sample collected at the T(4S') resonance with the Belle detector at the KEKB asymmetric- 
energy e+e~ collider. We fully reconstruct one neutral B meson in the p^7r^ final state and identify 
the fiavor of the accompanying B meson from its decay products. We obtain the charge asymmetry 
A'^p = — 0.16ib0.10(stat)ib0.02(syst). An unbinned maximum likelihood fit to the proper-time 
distributions yields Cp^ = 0.25±0.17(stat)l°;°^(syst), AC^^ = 0.38±0.18(stat)l°;°^(syst), Sf„ = 
-0.28±0.23(stat)+2:o8(syst), and ASp^, = -0.30±0.24(stat)±0.09(syst). The direct CP violation 
parameters for B —>■ p+7r~ and B —>■ p~iT'^ decays are = — 0.02±0.16(stat)lQQ2(syst) and 

Ap+ = -0.53±0.29(stat)tg:°^(syst). 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Gv, 14.40.Nd 
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In the Standard Model (SM) of elementary particles, CP violation arises from the 
Kobayashi-Maskawa (KM) phase in the weak-interaction quark-mixing matrix. Recently, 
the Belle i2| and BaBar collaborations [3] reported results on CP violation via h uud tran- 
sitions m. B TT^TT^ decays, which are related to the CP violation parameter 02- Here we 
present a study of 5 — > p^7r^, which is another 02 related decay. Since B p^vr^ is not 
a. CP eigenstate decay, four decay modes with different charge and flavor combinations in 
the neutral B system must be considered. 

In the decay chain T(45') — > — > (p^vr^)/tag, one of the B mesons decays at time 
tp-K to p'^n^ and the other meson decays at time ttag to a final state /tag that distinguishes 
between B^ and B^. The decay rate for B^{B^) —>■ p^Tr^ has a time dependence given by 

vt^^m = (i±^s>)^- — (1) 

x{l + q-[{Sp^ ± ASp^) sin(Am^At) 
-{Cp^ ± ACp^) cos(AmdAt)]}, 

where r^o is the B^ lifetime, Am^ is the mass difference between the two 5° mass eigenstates. 
At = tpTT — ttag, and the 6-fiavor charge q = 1) when the tagging B meson is a B^{B^). 
The time and flavor integrated charge asymmetry A'^(^p is defined as 

_ N{p+iT-) - N{p-7i+) 

N{p+n-) + N{p-7r+y ^ ^ 

where N{p~^7T~) and N{p~7r~^) are the sum of the yields for B^ and 5° decays to p~^7r~ 
and p~vr+, respectively. The mixing- induced CP violation parameter S'p^ is related to 02 
and CpTT is the flavor-dependent direct CP violation parameter. The asymmetry between 
the decay rates, 1(5° ^ p+vr-) + T{W p-n+) and 1(5° ^ p-n+) + T{W P+tt"), 
is described by ACp^r, while the strong phase difference between the amplitudes contribut- 
ing to B^ — > pvr decays is described by AS'p^. These parameters are related to 02 as 
Sp^ ± ASp^ = ^l-{Cp^±ACp^ysm{2cPis ± 5), where 20±g = arg[(g/p)(A± /AJJ] and 
S = arg[y4~^/A+^]; arg[g/p] is the B^-B^ mixing phase. The terms ^+^(^4+^) and A~^{A~^) 
denote the transition amplitudes for the processes B^{B^) —>■ p+7r~ and B^{B^) P^tt^, 
respectively. The angles 0^^^ are equal to 02 if there is no penguin contribution. The effect 
of direct CP violation can also be expressed in terms of another set of parameters, A'l^~ and 

iV(50 ^ pTvr±) + N{BO p±7rT) ^ 

_^-^CP ^ p-K ^ -^CP '-^^ p-K 

1 ih AC p-^ ih A^Qp • C p-j^ 

The strategy of this analysis is to reconstruct final states in quasi-two-body decays B^ — »• 
^^i^o^^T^ which correspond to distinct bands in the tt+tt^tt^ Dalitz plot. We exclude the 
interference region where the p charge is ambiguous, and neglect possible residual interference 
effects. 

The results for this analysis are based on 140 fb~^ of integrated luminosity, which corre- 
sponds to 152 X 10^ produced BB pairs. The data were collected with the Belle detector 
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at the KEKB asymmetric-energy e~^e~ collider [5!], which collides 8.0 GeV e~ and 3.5 GeV 
e"*" beams. The T(4S') is produced with a Lorentz boost of /?7 = 0.425 nearly along the 
electron beamline. Since the and mesons are approximately at rest in the T{4S) 
center-of-mass system (CM), At can be determined from Az, the displacement in z between 
the p^TT^ and /tag decay vertices: At ~ (-2p±7rT — Ztag)/P'yc. The z axis is anti-parallel to 
the positron beam. 

The Belle detector P] is a large-solid-angle general purpose spectrometer that consists of 
a silicon vertex detector (SVD), a central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), time-of- flight scintillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux return located outside of the coil is 
instrumented to detect mesons and identify muons. 

To reconstruct B^ —>■ p^vr^ candidates, we combine pairs of oppositely charged tracks 
with tt" candidates. Each charged track is required to have transverse momenta greater than 
100 MeV/c in the laboratory frame. Charged tracks are identified as pions by combining 
information from the ACC, CDC and TOF. Electron-like tracks are rejected. The 7 energies 
for 7r° candidates are required to be greater than 50 MeV if the photon is detected in the 
barrel ECL (32° < 9 < 129°); otherwise, the energy is required to be larger than 100 MeV, 
where 9 denotes the polar angle with respect to the z-axis. The 7r° candidates are selected 
from 77 pairs with invariant masses in the range 0.118 GeV/c^ < My^ < 0.150 GeV/c^, 
and momentum larger than 200 MeV/c in the laboratory frame. In addition, we require 
I cos6'Jgj.| < 0.95, where 6*5°^. is defined as the angle between the photon flight direction and 
the boost direction from the laboratory system in the 7r° rest frame, and we require the 
of the 7r° mass-constrained fit to be less than 50. 

B meson candidates are reconstructed using the beam-energy constrained mass Mbc = 
a/ -E^eam ^ ^"^^ energy difference AE = Eb — -E'beam- The variables Eb and Pb 
are the reconstructed energy and momentum of the B candidate in the CM frame, and 
-E'beam IS the CM beam energy. The B candidates in the region with Mbc > 5.2 GeV/c^ and 
-0.3 GeV < AE < 0.2 GeV are selected. The signal region is defined as Mbc > 5.27 GeV/c^ 
and —0.10 GeV < AE < 0.08 GeV. The B —>■ p^vr^ candidates are formed from 3-body 
B TT^TT^Tr" decays with a vr^'^vr^ invariant mass in the range 0.57 GeV/c^ < Mn-±7rO < 0.97 
GeV/c^ and p helicity | cos9^^i\ > 0.5, where 9^^^ is defined as the angle between the charged 
pion direction and the B^ direction in the p rest frame. To avoid the region where the 
p'^7T~ and p~7r~^ contributions interfere, we exclude candidates with both M^+^o and M^— ^0 
smaller than 1.22 GeV/c^. Candidates with Mt^+tt- < 0.97 GeV/c^ are removed to avoid the 
region where the p'^n^ or p^vr^ bands overlap with p^ir^. 

To suppress the dominant e~^e~ qq continuum background (g = u, d, s, c), we form the 
likelihood ratio TZ = Cg/ ('^s+'^bkg), where Cg and £bkg are likelihood functions for signal and 
continuum events, respectively. We use a Fisher discriminant based on five modified Fox- 
Wolfram moments 0], and the CM flight direction of the B {9b) with respect to the 2;-axis to 
form the likelihood function. The signal likelihood Cs is determined from a GEANT-based 
Monte Carlo (MC), and £bkg is based on Mbc sideband data, Mbc < 5.26 GeV/c^. The 
continuum background is reduced by requiring TZ to be greater than 0.8. If there is more 
than one candidate in an event, we select the candidate with the smallest sum of the for 
the TT+TT" vertex fit and the vr" mass-constrained fit. 

The flavor of the accompanying B meson is identified from the decay products not as- 
sociated with the reconstructed B^ —>■ p'^Tc^ decay. We use the same method as used for 
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FIG. 1: AE (top) and M^c (bottom) projections for the result of the 2-D unbinned hkehhood fit. 
The plots on the left are the results for the p+vr^ candidates, while those on the right show the 
results for the p~Tr^ candidates. 



the Belle sin20i measurement [8|, |9|. Two parameters q and r are used to describe the 
flavor tagging information. The parameter q is defined in Eq. ^ and the parameter r is a 
MC-determined quality factor that ranges from r = for no flavor discrimination to r = 1 
for unambiguous flavor assignment. It is used only to sort data into six r intervals. The 
wrong tag fractions for and are obtained from B D*lu, D*n, D*p, and Dtt data 
for the six r intervals. 

The vertex reconstruction algorithm is the same as that used for the sin 20i analysis . 
The vertex positions for p^vr^ and /tag decays are reconstructed from charged tracks with 
associated SVD hits and an interaction point constraint. The vertex for /tag is determined 
from all well-reconstructed tracks excluding the tracks from the B^ —>■ p^Tr^ decay and Kg 
candidates. 

Figure Q shows the AE (Mbc) distribution in the Mbc {AE) signal region for B^ 
p^TT^ candidates after flavor tagging and vertex reconstruction. The p^7r^ signal yields 
are extracted from an unbinned maximum-likelihood fit to the two-dimensional (MbcAi?) 
distribution. The backgrounds are categorized as continuum qq, b —>■ c transitions [BB), 
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FIG. 2: Signal yields as functions of (left) Mtttt and (right) cos 9^^^ in data. The histograms show 
the results of i? ^ p^vr^ MC simulation with areas normalized to the total signal yield. 



B — > pK, and rare charmless decays other than B pK (rare B). The distributions for 
pvr, BB, pK, and rare B events are obtained from MC. 

The pn signal PDF contains two components: signal events reconstructed with the correct 
charge (Ppn) and those with incorrect charge {PJ^)- The fraction of events with incorrect 
charge in the signal region due to combinations that include a random 7r° is estimated to 
be 2.7% from MC and is fixed in the fit. The signal PDF shape is modeled by a smoothed 
histogram. The AE distributions for B ^ pn signal are parameterized separately for 7r° 
momentum below and above 1.2 GeV/c in the laboratory frame. The AE widths for pvr 
and pK are calibrated from D°[K-n+]n° data. The B+ D'^[K-n+n'^]7i+ mode 

is used to calibrate the AE and Mbc peak positions. The Mbc and AE distributions for 
the continuum qq are parametrized by an ARGUS background function ^] and a linear 
function, respectively. The contributions from B pK (with B = (9.0 ± 1.6) x 10~^ [lll |) 
and from rare B decays are fixed in the fit, while the yields for B ^ pn signal, BB and 
continuum backgrounds, and the shape parameters for continuum are floated. We obtain 
483 ±46 B 



p^ir^ events, and obtain a time and flavor integrated charge asymmetry 



pvr, B pK, qq, BB and rare 



—0.16 ± O.lO(stat). The estimated yields for B - 
B in the signal region are 328.7, 11.2, 833.0, 23.3 and 18.8, respectively. We remove the 
requirements on M^±^o and cos 9^^^ and examine these distributions to verify that the signals 
reconstructed as B ^ n~^7r~n'^ are from the two-body decay B pvr. Figure El shows the 
signal yields in bins of M^±^o and cos 9^^^ for data. 

The CP violation parameters are obtained from an unbinned maximum-likelihood fit to 
the observed proper-time distribution for the B ^ pn candidates in the (Mbc,Ai?) signal 
region. The likelihood function describing the proper-time distribution is 



N 



^ = l[{fp.PpA^u) + f;:p;:{Au 



(4) 



i=l 



+ fpKPpK{Ati) + fqqPqq{Ati) 

(At.)}, 

where the weighting functions fm {m = pn, pK, qq, BB, and rare B) are determined on 
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an event-by-event basis as functions of Mbc and AE for each flavor tagging r interval and 
tt" momentum range in the laboratory system. The time-dependent probability density 
functions {At PDFs) PpT,{Ati) for B ^ pn and Ppx(Atj) for B — pK are obtained from the 
true PDFs convolved with the At resolution function used in the sin20i measurement 0. 
The true PDF for B ^ pir is given by Eq.l modified to incorporate the effect of incorrect 
flavor tagging. The PDF for B ^ pir signal reconstructed with incorrect charge, P^^:p(Atj), 

is given by Pp^^±{Ati). Foi B ^ pK, C = S = AS = 0, AC = -1, and A^^ = is 
assumed. The resolution function consists of the detector resolution, the shift in vertex 
position due to secondary tracks originating from charmed particle decays, and smearing 
due to the approximation At ~ (-2p±7rT — Ztag) / Pjc. The At PDFs for other backgrounds are 

all parameterized as Vj = {1 — fj)6{At — p^) + /-,exp(— ) convolved with Rj {j = qq, 

BB and rare B), where fj is the fraction of the background with effective lifetime tj. The 
resolution-like function Rj for background is given by two Gaussians. The parameters of the 
At PDF for qq background are obtained from a fit to sideband data (5.2 GeV/c^ < Mbc <5.26 
GeV/c^ and AE > —0.15 GeV). The parameters of the At PDFs for BB and rare B are 
obtained from a fit to MC. 

The maximum likelihood fit to the 1,215 pir candidates gives C^tt = 0.25 ± 0.17l^oQg, 
ACp^ = 0.38 ± O.lStoM, Sp^ = -0.28 ± 0.23+[J:J^ and ASp^ = -0.30 ± 0.24 ± 0.09, where 
the first (second) errors are statistical (systematic). The correlation between Cp-,, and ACp^ 
is 0.271 and that between Sp^ and ASp-jr is 0.284, while correlations between other variables 
are smaller. The data and fit result are shown in Fig. El 

The systematic error in A'^p includes a possible background asymmetry (±0.010) and 
charge asymmetry in the tracking (±0.012). The charge dependence of tracking efficiency 
is studied using — > K^n~^ decays from inclusive D*~^ D^7r~^ and selecting the momen- 
tum region corresponding to B^ p'^n^ decays. The systematic errors for time-dependent 
measurements include the uncertainties in the vertex reconstruction, background fraction, 
background At PDF, wrong-tag fractions, pir and pK At resolution functions, physics pa- 
rameters {tb, Arrid 3|; A''^ ^^) and fitting bias. The fitting bias is estimated from MC 



pseudo-experiments. All other systematic uncertainties are obtained by varying parameters 
within their errors and repeating the fit. The dominant source of systematic error is the 
vertex reconstruction for C^^, 1°:°^^ for ACp^, 1°:°?^ for Sp^, and for ASp^). 

We perform various consistency checks. We examine the stability of the results as the 
TZ selection criterion is varied and the asymmetry of the At distributions for events in the 
sideband region. No significant variation or asymmetry is observed. We measure the B^ 
lifetime with the B^ — > p^n^ candidates and find r^o = 1.56lg'i2 PS? which is consistent 
with the world average value 1^ |. 



The extraction of 02 from measurements of time-dependent CP violation parameters 
in i? ^ P^tt'f decays has been studied in several theoretical approaches 3, Q Q. A 



Grossman- Quinn type bound 16] based on isospin (SU(2) symmetry) does not signifi- 
cantly limit the penguin diagram contribution due to the large branching fraction for 
B^ — > p^TT*^ 17]. Since the number of measurable quantities (six including B{B^ —>■ p~^7r~)) 
are not sufficient to completely describe the amplitudes for B^ —>■ p^ir"^ decay (8 free param- 
eters), either specific models or additional assumptions are involved, such as QCD factor- 
ization or SU(3) fiavor symmetry |^. A recent approach assuming broken fiavor-SU(3) 
implies 02 = (102 ± 11 ± 15)° using our results The first error is experimental while 
the second is the uncertainty due to SU(3) breaking effects. 

In summary, using 152 x 10^ BB pairs, we have measured CP violation parameters 
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FIG. 3: At distributions for —>■ p^vr^. (Top) B° and tagged p'^TT and p 7r+ candidates. 
(Bottom) raw CP asymmetries in high and low r intervals for p+vr^ and p^ir^ . The solid curves 
show the fit results. 

for fiO ^ p^n^ decays. We obtain A^Jp = -0.16 ± 0.10 ± 0.02, Cp^ = 0.25 ± 0.17;|J-|J^, 
ACp^ = 0.38±0.18^[j:[J^, Sp^ = -0.28±0.23lS;J^ and ASp^ = -0.30±0.24±0.09. These give 
the direct CP violation parameters A^' = -0.02 ±0.16l^^ and = -0.53 ±0.29^°;°^. 
These results are consistent with a previous measurement [13|. We find no significant mixing- 
induced or direct CP violation in — > p^vr^. 
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